1. Joint Modeling
1.1 Trace plots for convergence check
The current MCMC setting is:
- 100,000 iteration;
- 90,000 burn-in;
- 10 thinning.
1.2 Gelman and Rubin’s convergence check
## Potential scale reduction factors:
##
## Point est. Upper C.I.
## HDevsum 1.01 1.04
## LDevsum 1.01 1.02
## dh0 1.42 2.05
## dh1 1.13 1.36
## dl0 1.51 2.22
## dl1 1.06 1.16
##
## Multivariate psrf
##
## 1.61
1.3 ACF Plots
Here we plotted ACF plots for the following variables:
- Total deviance;
- Variables that didn’t pass the convergence check.
1.4 WAIC results
| LevelH | LevelL | |
|---|---|---|
| DIC | 1343.40738 | 22794.8154 |
| DIC3 | 1218.27290 | 22823.8719 |
| PWAIC | 63.37221 | 196.0088 |
| WAIC | 1258.22345 | 22835.7281 |
2. Separate Modeling
2.1 Trace plots for convergence check
The current MCMC setting is:
- 100,000 iteration;
- 90,000 burn-in;
- 10 thinning.
2.2 Gelman and Rubin’s convergence check
## Potential scale reduction factors:
##
## Point est. Upper C.I.
## HDevsum 1.00 1.00
## LDevsum 1.00 1.00
## dh0 1.84 3.05
## dh1 1.86 2.98
## dl0 1.28 1.71
## dl1 1.09 1.26
##
## Multivariate psrf
##
## 1.71
2.3 ACF Plots
Here we plotted ACF plots for the following variables:
- Total deviance;
- Variables that didn’t pass the convergence check.
2.4 WAIC results
| LevelH | LevelL | |
|---|---|---|
| DIC | 1353.59804 | 33770.3243 |
| DIC3 | 1230.04462 | 33889.3234 |
| PWAIC | 84.87555 | 167.7255 |
| WAIC | 1295.26122 | 33894.1349 |